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ABSTRACT 



We present a method using the SALT2 light curve fitter to determine the redshift of Type la supernovae in the 
Supernova Legacy Survey (SNLS) based on their photometry in g' , r' , i' and z . On 289 supernovae of the first three 
years of SNLS data, we obtain a precision <Jaz/(i+z) = 0.022 on average up to a redshift of 1.0, with a higher precision 
^ ' of 0.016 for z < 0.45 and a lower one of 0.025 for z > 0.45. The rate of events with |Az|/(l + z) > 0.15 (catastrophic 

Q\ ' errors) is 1.4%, and reduces to 0.4% when restricting the test sample of spectroscopically confirmed Type la supernovae. 

Both the precision and the rate of catastrophic errors are better than what can be currently obtained using host galaxy 
• photometric redshifts. Photometric redshifts of this precision may be useful for future experiments which aim to discover 

up to millions of supernovae la but without spectroscopy for most of them. 



Key words. Supernovae 

1. Introduction indeed have to count on an alternative to select the SNe la 

and determine their redshift from light curves only. 



Ten years after the discovery of the accelerating expan- Studies solely based on photometry present the advan- 

H \ sion of the Universe, Type la supernovae (SNe la) are still tage of leading to larger sets of supernovae, that can be 

ri , among the most accurate probes of dark energy. To study use d to compare properties of SNe la in various subsam- 

systematic effects more precisely, future experiments aim p l es , sorted for instance by redshift, SN color, or host mor- 

at an increase of several orders of magnitudes in the num- phology. The possibility to use solely photometric redshifts 

ber of detected supernovae. Pan-STARRS Medium Deep to derive cosmological constraints from SNe la, however, 

Survey (MDS) and the Dark Energy Survey will discover remains to be shown, 
thousands of supernovae up to a redshift of 1, though they 

will only be able to observe spectroscopically a fraction of This paper addresses the question of obtaining photo- 

them. The Large Synoptic Survey Telescope will discover metric redshifts for SNe la. It presents a method that can 

roughly 250,000 Type la SNe per year in the same redshift be applied to ongoing as well as to future projects. The 

range, but again with o nly partial sp e ctrosc opic coverage combination of a photometric selection of supernovae in 

(see the white paper of lHowell et all (120091) for a review the SNLS data, which doubles the number of supernovae 

of future survey s). While current su rveys as the Supernova compared to the real-time selection, with the photometric 

Legacy Survey (|Astier et all I2006D are already limited by redshifts derived according to the method presented in this 

the time available for spectroscopy, most future surveys will paper, will be the subject of a forthcoming study. 



2 



N. Palanque-Delabrouille et al.: Photometric redshifts for SNIa in SNLS 



Previous works have already tackled the topic of SN la 
photometric redshifts. An empirical p hotometric redshift 
estimate for SNe la was proposed by IWang et al.l (|2007|) 
who obtained a rms dispersion of 0.050 for 20 supernovae 
from the SNLS first year data not contained in the train- 
ing set, and 0.031 for 20 supernovae used in the training. 
A method was also developed in SNLS for the needs of the 
real-time selection of SN la c andidates to send for spec- 
troscopy (|Sullivan et all [2006). It used the rising part of 
the light curves and a cosmology prior to both select ten- 
tative SN la candidates and derive an estimate of their 
redshift. Fits to the pre-maximum light data reached a pre- 
cision a\Az\ /(t+z) — 0.046 while fits to the entire light curve 
a precision of 0.037. 

The method presented in this paper aims at an im- 
proved precision on the photometric redshift, while remain- 
ing independent of the underlying cosmology. The study 
was calibrated with simulated SN la light curves, whose 
characteristics are given in §2. The method is presented 
in §3. To test its performance, it was applied to a large 
sample of SNe la from the first three years of SNLS data, 
presented in §2. The results obtained on this test sample 
are given in §4. A discussion on the method and on the 
impact of the underlying cosmology is presented in §5. The 
possible redsfhift biases and the precision of the method for 
various sub-classes of SNe la is discussed in §6. Conclusions 
are then given in §7. 

2. SNLS data and simulation description 

The SNLS data cover the redshift range between and 
1.2. They consist of four-band light curves obtained in a 
rolling search mode, leading to roughly 30 points per ob- 
serving season in r' and i', 20 points in g' and 15 points 
in z' . The test sample used to quantify the perfomance of 
the method consists of a total of 289 SNe la passing color 
and stretch cuts (|c| < 0.4, \x\\ < 4) to avoid events with 
extreme values compared to typical SNe la. It originates 
from the union of two sources: 241 are real-time selected 
events confirmed as "la" o r "la?" by spectroscopy (|Howel 
2005; iBronder et all 120081: lEllis et all . 120081: iBalland et a" 



ure[T|). These distributions best reproduce our selected data 
sample at redshift z < 0.7 where Malmquist bias is negligi- 
ble (cf.Hi^ir^] (HOOK) fig. 13). 



2009), and the additional 48 are issued from a SN la pho- 
tometric selection by Bazin et al. (in prep.) based on the 
host galaxy photometric redshift. The latter SNe la have a 
host spectroscopic redshift (thus allowing their use in the 
test sample) but no type determination because of insuffi- 
cient S/N for the supernova spectrum. The contamination 
of the photometric selection of SNe la is estimated to be 
7%, implying a possible contamination 3.4 events out of 
48. Our total sample of 289 events is thus expected to be 
contaminated at the 1% level. 

Synthet ic SN l a light curves were produced with SALT2 
(|Ouv et all 12007ft assuming a flat Universe with Qm — 
0.26. The light curves are characterized by five parameters: 
color c, stretch Xi, absolute normalization, date of maxi- 
mum flux in the rest frame B band and redshift. Redshifts 
were simulated in the range from to 1.2 assuming a con- 
stant volumetric SN la rate. The date of maximum was 
drawn uniformly within the dates of observation in SNLS. 
Values of stretch were generated according to a Gaussian 
distribution (centered on m = and with standard de- 
viation (T = l), and those of color according to a double 
Gaussian distribution (mi = 0.00, a\ = 0.05, m 2 — 0.08, 
02 = 0.12, and relative amplitudes A1/A2 = 0.91, see fig- 
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Fig. 1. Color distribution of SNe la with redshift z < 0.7 
from the first three years of SNLS data. The curve is a fit to 
the data with a double Gaussian as explained in the text. 



Model light curves in the four bands of SNLS (</, r' , i' 
and z') were computed. Instrumental effects were taken into 
account following the distributions of flux uncertainties and 
dispersions observed in data, as well as correlations between 
flux uncertainties and flux values. A total of 20000 SNe la 
were generated. A preselection filter was applied to these 
events to remove those where the light curve was either too 
faint or too poorly sampled to be distinguished from pure 
noise, leaving 11106 SNe la. More details on the simulation 
and the preselection can be found in Bazin et al. (in prep.). 



3. Redshift fit 

The light curve fitter used for this analysis is based on the 
public version of SALT2 with additions described below. 
The x 2 minimization is done with the Minuit packageQ 

By default, the fitting procedure deals with redshift as 
with the other free parameters, using initial values of 0, 
and ~ 0.1 respectively for the parameters describing the 
color c, the stretch x\ and the redshift of the supernova. 
The fit consists of a x 2 minimization of the multi-band 
light-curves according to an empirical modelling of the SN 
la luminosity variations as a function of phase, wavelength, 
X\ , and c. Figure [5] shows the distribution of photometric 
redshift (z p ho) vs. spectroscopic redshift (z spe ) obtained us- 
ing the default SALT2 to fit z p ho- Despite the use of priors 
on color and stretch defined by the statistical distribution 
of these two parameters (see section l3~TT) . 43% of the events 
had |z p ho-.2spc|/(l + Zspo) = |Az|/(l + z) > 0.07, and some 
redshifts, such as z p h ~ 0.4, were clearly favored over oth- 
ers. This shows that an initial step is required in order to 
start the x 2 minimization with reasonable values of the pa- 
rameters for each individual supernova (not surprisingly, as 
the fit is non- linear). The procedure to initialize the fit is 
described in section [ 



http:/ /wwwasdoc. web. cern.ch/wwwasdoc/minuit/mmmain. html 
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Fig. 2. Photometric vs. spectroscopic redshift for the data 
sample in the absence of initialization or prior on redshift. 

The best precision one can achieve with SNLS data can 
be estimated by initializin g the fit at the a ctual value of 
the redshift, as was done in lGuv et al.l ((2007). This led to a 
precision (cf. definition in sectional) of 0.018 and only 3% of 
events with |A|z/(l + z) > 0.07. The main goal of the study 
presented in this paper was to approach this ideal situation 
by determining, prior to the full fit, the most accurate initial 
values not only for redshift but also for color and stretch. 

Since different combinations of the five parameters that 
describe the four-band light curves of a SN la could yield 
similar x 2 (although giving extreme values of color or 
stretch) , it was essential to guide the fitter with priors in or- 
der to converge to the most accurate values (cf. section l3~Tj) . 

3.1. Priors and bounds 

Several of the parameters pj were constrained by Gaussian 
priors (mean rrij and variance Vj ) which were included in a 
generalized x 2 as follows: 

[/data(^z) /model(^i)] 



E 



E 



(p.1 



(1) 



where f(U) is the flux at date and cr(tj) is the flux un- 
certainty at the same date. 

SNe la form a fairly homogeneous class of objects. 
Significant departures from typical values of color and 
stretch could therefore be avoided by setting Gaussian pri- 
ors on these two parameters: mean m c — 0.04, standard 
deviation er c = 0.15 for color and m Xl = 0, cr Xl = 1 for 
stretch, in agreement with the data. 

Redshift bounds and prior were estimated in the fol- 
lowing way. The four band light curves were fitted by the 
phenomenological form 

p —(t— to)/TfaU 

/(*) - V e -( t -,o)/^ +B ^ 

to obtain estimates of the observed peak luminosities. As 
SNe la are standard candles, a strong correlation is ob- 
served between the redshift and the estimated peak flux 
/max determined from the fit by equation [2J This correla- 
tion, illustrated in figure [3] for the r' band, allowed setting 
lower (zi) and upper (22) bounds on redshift as well as es- 
timating a mean expected redshift value m z , along with an 



uncertainty <j z taken as the r.m.s. dispersion around m z for 
a given observed / ma x- Typical values for o z lie between 0.06 
and 0.11. Two SNe la lie outside the limits shown in the 
figure. One has poor sampling (the first point after peak in 
any of the four bands occurs after 25 days). The other has 
both a low stretch (x\ — —1.2) and a large intrinsic color 
(c = 0.36, close to our limit of 0.4 for our data sample). 




4.5 5 
log (f max) 

Fig. 3. Correlation of redshift with the peak flux in r' . 
Black points are for simulated SNe la, red dots for the 
data sample. The blue curves represent the lower and up- 
per bounds set on redshift, and the green curve the mean 
estimate of redshift for a given peak flux. 



Independent estimates of the redshift bounds z\ and 
Z2j mean value m z and dispersion a z were computed for 
the two bands with the best sampling and photometric ac- 
curacy, namely r' and i'. As a prior on redshift, we used 
a Gaussian with mean set to the weighted average of the 
m z values from r' and i' , and with variance taken to be 
[o\z(i') -2 + cr 2 (r') -2 ] -1 . The bounds being quite conserva- 
tive (over 99.5% of the events contained within the limits), 
the most restrictive ones were retained to limit the range of 
redshifts probed. The bounds are only useful to reduce the 
time consumption of the fit and do not lead to further im- 
provement beyond what is already achieved with the prior 
on redshift. 



3.2. Fitting procedure 

In the standard use of SALT2, the fit of the light curves is 
done at known redshift, and it is best to use only the bands 
which fall in the rest-frame wavelength range between 290 
and 700 nm. In the case studied here, as the redshift is a 
free parameter, all four bands g', r' , i 1 and z' are used at 
all times. The priors on color, stretch and redshift defined 
above were used throughout. 

The first step of the fit consisted in computing the re- 
duced x 2 f° r successive values of the redshift separated by 
Sz = 0.1, within the bounds defined as explained above. 
To reduce the time consumption, this step is done at fixed 
c = m c and X\ — m Xl . The value of the redshift z m in|c,:ri 
yielding the minimum reduced x 2 m this scan was taken 
as the initial redshift for the second step. In the presence 
of two minima, both were considered thereafter. Note that 
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with the bounds set on the redshift range, multiple minima 
only occurred in 3.5% of the events. 

A finer tuning of the initial redshift was then obtained 
by setting the color and stretch parameters free and per- 
forming a new redshift scan around z m i n |c,a:ij with the same 
step 8z as above, until a reduced x 2 minimum was reached. 
The values of redshift, color and stretch for the fit that 
yielded the lowest reduced \ 2 were taken as initial values 
for the full fit. The latter was performed with the stan- 
dard code of the light curve fitter SALT2, letting all five 
parameters free and initialized as described above. When 
two minima were observed in step 1, both were taken into 
account. The solution with the smallest reduced \ 2 at the 
end of the full fit was chosen. Note that in 71% of the cases 
with two x 2 local minima, both minima ended up with the 
same set of values for all five parameters (within statisti- 
cal errors). Only 1% of the events thus truly needed to go 
through all 3 steps for both z m i n | C 2;i to determine the best 
redshift. 



4. Statistical accuracy of SN la photometric 
redshifts 

The performance of the method described above was as- 
sessed from the two samples of simulated and observed 
SN la light curves. The photometric redshift was compared 
with the generated one (z gGn ) for simulated events and with 
the spectroscopic one (z spe ) for data. 

In both cases, 99% of the events lie within 0.15 x 
(l + z) of the actual redshift. With the precision de- 
fined as in ICoupon et al.1 (|2009T ) by <Jaz/(i+z) = 1-48 x 
median[| Az|/(1 + z)] (a robust estimate of the r.m.s), and 
the outlier rate 77, or rate of catastrophic errors, defined as 
the proportion of events with |Az|/(l + z) > 0.15, we ob- 
tain a precision of 0.022 and 77 = 1.4% on average over the 
entire redshift rang^3 (see figure I3J. 10% of the events have 
|Az|/(l + z) > 0.07. For 70% of these, the reduced x 2 is 
smaller than for the fit initialized at the actual redshift^ 
The photometric redshift therefore appears, wrongly, to be 
a better solution than the actual redshift. This is a limita- 
tion of any method based on \ 2 minimization. 

The accuracy achieved for two redshift intervals (below 
and above 0.45) is given in table [T]. The reduced accuracy 
at high redshift comes from the vanishing SN flux in the 
g' then r' bands of SNLS, thus enhancing fit degeneracies 
(see discussion) . The redshift dependence of the accuracy is 
well reproduced by an increase of the spread in Az propor- 
tionally to (1 + z) 2 , with <jaz/(i+z) 2 = 0.014 independently 
of redshift . 

Comparison of columns 2 and 6 shows that 77 for sim- 
ulation and data are in agreement, given the uncertainty 
of 0.7% on 77 for data. Resolutions at high redshift are also 
in agreement. At low redshift, it is slightly worse for data, 
probably due to sources of SN la heterogeneity not repro- 
duced in the simulation. 



2 SALT2 was trained with 40% of the data used here. 
Restricting to the complementary set yields results in agreement 
with the quoted performances within statistical uncertainties. 

3 Note that the \ 2 comparison is done on the sole contribution 
of the data points (i.e. not including the contribution of the 
priors), as there is no redshift prior for the fit initialized at the 
actual redshift. 
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Fig. 4. Photometric vs. spectroscopic redshift for the 
method of this paper. The green lines are for z p h = 
z sp o ± 0.022(1 + z spo ) representing the average precision 
over the entire redshift range, and the blue lines for Zph G — 
z sp o i 0.15(1 + z spe ) to visualize the catastrophic redshifts. 
Red circles are data SNe la for which a spectroscopic red- 
shift is available, black points are the simulated SNe la. 

As visible in columns 5 and 6, the photometric de- 
termination of the supernova redshift has very few out- 
liers for spectroscopically confirmed SNe la. Most of the 
catastrophic failures occur for events that have been added 
by the photometric selection. The increase in 77 between 
columns 5 and 6 is in agreement with the estimated con- 
tamination fraction of 1% of the total sample. The preci- 
sion also tends to be slightly better for spectroscopically 
confirmed events. 

An alternative to supernovae photometric redshifts is 
to use pho t ometr ic redshifts of the host galaxies, z ga i. 
Illbert et al.l (|2006f ) have published photometric redshifts for 
522286 objects in the four deep fields of SNLS. Only 83% of 
potential SNe la have a matched hos10 with a z ga j, while SN 
la photometric redshifts can be computed for all objects. 
The two photometric redshift determinations are compared 
in the last two columns of tabic [TJ Supernova photometric 
redshifts have a much better precision. However, when com- 
bined with a photometric selection of SNe la, it may lead 
to samples more contaminated by core-collapse supernovae 
for instance than when using galactic photometric redshifts 
which are independent of the light curves. Future projects 
might want to combine the two techniques to obtain a se- 
lection of SNe la with reduced contamination from non-la 
supernovae (using fits with redshift fixed at z ga i to obtain 
the selected sample of SNe la) while having the redshift 
accuracy of supernovae photometric redshifts (refitting the 
redshift of the selected events to derive z p ho)- 

5. Discussion 

As illustrated in columns 1 and 2 (for the simulation) or 
4 and 6 (for the data) of table [TJ the prior set on redshift 

4 The match to a host galaxy was considered successful out to 
a supernova- host distance of 5 r ga i with r ga i defined as the half- 
width of the galaxy in the direction of the event. 9% of the SNe 
lie in regions masked in the galaxy catalog because of neighbor- 
ing bright stars, and 9% of the remaining SNe are associated to 
a galaxy lacking a secure z ga i for diverse photometric reasons. 
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(2) 


(3) 


(4) 


(5) 
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(7) 


ffA*/(i+«) for z 


< 0.45 


0.006 


0.006 


0.018 


0.019 


0.014 


0.016 


0.031 


"a 2 /(i+ 2 ) for z 


> 0.45 


0.029 


0.027 


0.024 


0.027 


0.020 


0.025 


0.037 


Outlier rate 77 




2.5% 


1.0% 


2.0% 


3.1% 


0.4% 


1.4% 


5.5% 



Table 1. Accuracy 0^2/(1+2) of the supernova photometric redshift z p h : columns 1 and 2 are for simulated events, 
columns 3 to 6 for data. For comparison, column 7 indicates the accuracy of the photometric redshift of the host galaxy 
Zg a i. The data sample contains 57 objects at z spc < 0.45 of which 50 are spectroscopically confirmed, and 232 at z spc > 0.45 
including 191 spectroscopically confirmed SNe la. Superscript b if bounds and prior on redshift used (as explained in 
section [3]), a otherwise. 



reduces the rate of catastrophic errors by ~ 1.5%. Their 
impact on precision, however, is small, because the precision 
attained is significantly better than the width a z of the 
prior. 

We have also checked the impact of assuming that SNe 
la are "normal" by removing all the priors on color, stretch 
and redshift but maintaining the same fitting procedure, 
i.e. two consecutive \ 2 scans to initialize the values of the 
parameters for the full fit. The precision obtained on aver- 
age over the entire redshift range is then of 0.029 and the 
rate of outliers is significantly increased to 77 = 9.6%. Most 
of these accidental fits end with extreme colors (c > 2) as- 
signed to the events. The assumption that SNe la are nor- 
mal therefore improves slightly the accuracy (by ~ 25%) 
but is crucial to reach a small rate of accidental redshifts 
(reduction by factor ~ 10). 

As inaccurate redshifts were correlated with inaccurate 
colors, we tried to replace the statistical prior on color by 
an estimate of the color of each event based on its four-band 
photometry. This led to little improvement due to the large 
dispersion around this estimate (often larger than a c = 0.15 
observed statistically on the sample of data SNe la). We 
thus remained with the prior m c = 0.04 and <r c = 0.15. 

Despite the initialization procedure through two consec- 
utive x 2 scans, some redshifts are still favored over others. 
To understand the origin of this problem, two tests were 
performed. 

The first one consisted in applying the procedure to per- 
fect simulated light curves with a sampling of one point per 
night and no flux dispersion. The redshift was determined 
with an accuracy <7a*/(i+«) < 0.001 over the full redshift 
range with no redshift favored (figure [5J top plot). 

In the second test, the light curves still had perfect pho- 
tometry but the sampling was degraded to reproduce that 
of real data. While the accuracy remained excellent, some 
accumulations (bottom plot of figure [5]) appeared at the 
same values of z p h as with the full simulation of figure 2] 
The z ~ 0.47 and z ~ 0.80 accumulations correspond to 
turn-overs in the r — i vs. redshift diagram (cf. figure [6]). 
The smaller one at z ~ 0.60 corresponds to a short plateau. 
Since SALT2 fits light curve colors, and as the r' and i' 
bands are the ones which contribute the most (best sam- 
pling and S/N), double redshift values for a given r — i color 
can indeed lead to the observed degeneracies. The first two 
can partially be broken with data in the g' band since the 
g — r vs. redshift diagram is monotonic (this also explains 
why the turn-over at z ~ 0.3 does not lead to redshift ac- 



cumulations). This is not the case for z ~ 0.8 where there 
is no longer any significant flux in the g' band. 

Note that when reaching perfect photometry, the data 
X 2 vanishes for the generated values of the fitted parame- 
ters and the priors then play a major role. In the case of 
redshift degeneracies, the fit converges to the redshift (out 
of the two possible ones) that minimizes the contribution of 
the color and stretch priors. Thus, near z p ho ~ 0.8, events 
with large/small intrinsic colors tend to be fitted with over- 
/under-estimated redshifts. With more data points in a 
third band (i.e. other than r' and the contribution of 
the priors to the total x 2 becomes negligible and the fit is 
truly constrained by the data only. As shown in the top 
plot of figure [5j a sampling of 1 point every day in all four 
bands fully breaks the degeneracies. 

As illustrated in figure [JJ the main impact of the pho- 
tometric uncertainties is to reduce the accuracy with which 
the redshift is determined by increasing the width of the 
Zpho ~ z spc band. Note that the distributions are non- 
Gaussian, which justifies the use of a median-based esti- 
mate of the r.m.s (section 0}. We have verified that the 
region thus defined does include at least 68% of the simu- 
lated events whatever the redshift range up to z ~ 1.1. 

SALT2 requires no assumption on cosmology. A depen- 
dence on cosmology could only come from the bounds and 
prior on redshift. This was tested by simulating SNe la in 
a flat fijif = 1 Universe where the SNe la are more lumi- 
nous by ~ 0.5 mag at z = 0.6. The bounds are sufficiently 
loose to be valid even in that extreme case. The prior on 
redshift is systematically underestimated by ~ 0.1, but at 
the end of the fitting procedure we obtain no bias and the 
same precision and outlier rate as in table ED column 2. The 
only difference is that outliers are more often under- than 
over-estimates of the redshift. 



6. Uncertainties and biases 

Although the photometric redshift shows negligible bias on 
average ((z p ho — z spc ) = 0.003), there is a clear correlation 
between color and redshift biases: overestimated redshifts 
are systematically correlated with underestimated colors. 

Supernovae with z < 0.6 have non-zero flux in at least 3 
bands and usually end up with good estimates of both red- 
shift and color. As the redshift increases however, the loss of 
information (due to unmeasurably small fluxes in the bluest 
bands) gives increasingly more weight to the color prior 
compared to the data. The fitted color thus tends to con- 
verge to that of the prior (m c ), leading to a color-dependent 
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Fig. 5. Photometric vs. generated redshift for simulated 
supernovae with perfect photometry. Time sampling is one 
point per day (top plot) or as in SNLS data (bottom plot, 
with a different scale). Lines indicate redshifts which tend 
to be favored by the fit. Results from sampling with one 
point per day superimposed on bottom plot as open green 
circles. 




Fig. 6. Generated colors as a function of generated red- 
shift. Shown in each plot are the events with a magnitude 
brighter than 25 in the two concerned filters. Lines indicate 
redshifts which tend to be favored by the fit. 



Fig. 7. Accuracy of the redshift determination in the full 
simulation (solid curve) and in the simulation with per- 
fect photometry but time sampling as in the data (dashed 
curve). Mean and rms are given for the full simulation. 

bias (cf figure [5]). This is particularly clear near redshift de- 
generacies (cf figure 0]), where events with positive (resp. 
negative) intrinsic colors populate the z p ho > z gcn (resp. 
Zpho < Zgcn) accumulations. 




-0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 

generated color 

Fig. 8. Redshift bias as a function of the supernova intrin- 
sic color for three redshift bins: red squares for z < 0.45, 
filled black circles for 0.45 < z < 0.65 and open blue circles 
for z > 0.65. The vertical lines delimit the la color range. 



To reduce this bias and improve the color and redshift 
determination for "colored" SNe (|c gcn — m c | > 0.1 where 
Cgen is the generated color) in SNLS, we performed an addi- 
tional SALT2 fit as follows : the mean and a of the Gaussian 
prior on x\ are set to the best fit value of x\ and its un- 
certainty from the procedure of section [31 a Gaussian prior 
on redshift is used with mean set to the value found previ- 
ously and er taken as 0.014 x (l + z p ho) 2 - With these stronger 
constraints, we eliminate the prior on color. 
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Above a redshift of 0.65, the additional fit reduces the 
bias Zpho — z gcn by ~ 0.02 for blue SNe (c gon — tti c < —0.1) 
while yielding a smaller but systematic improvement on 
red (cgen — rn c > 0.1) SNe. The color determination is also 
systematically improved: above a redshift of 0.65, the color 
bias is reduced by 0.04 for blue SNe and by 0.01 for red ones. 
Fits to events with c gen ~ ii c (black squares in figure [SJ , 
or to low redshift events where the data have more weight 
than the priors, are unchanged: these events were already 
well fitted. 

Figure IH1 illustrates the remaining bias on fitted redshift 
or fitted color as a function of redshift, for several ranges of 
intrinsic color. At high redshift, the reddest (and thus also 
the faintest) SNe exhibit the largest bias, emphasizing the 
lack of another filter (beyond the z band) with measurable 
flux to determine the SN parameters. On average for events 
passing the preselection filter, however, the bias remains 
small, with z p h — z gcn < 0.027 whatever the redshift up 
to a z r~j 1.1, and (z p h — z gcn ) = 0.008 on average over 
the full redshift range between and 1.1. The statistical 
uncertainty derived in section 2] is indicated in figure [§] It 
encompasses about 68% of the events at all redshifts. 




zpho 



accuracy reaches <J& z /(i+z) = 0.022 on average over the en- 
tire redshift range, with <7a z /(i+z) = 0.016 for z < 0.45 and 
°A*/(i+*) — 0.025 for 0.45 < z < 1.1. Degeneracies in the 
fit above z ~ 0.45, due to vanishing flux in the bluest band 
of SNLS, are responsible for the lower accuracy at high red- 
shifts. The spread in Az/(1 + z) 2 is seen to be independent 
of redshift, with oa 2 /(i+z) 2 = 0.014. 

Red (resp. blue) SNe tend to have an overestimated 
(resp. underestimated) redshift. The net bias however on 
the fitted redshift, after integration over the SN popula- 
tion, remains smaller than 0.027 whatever the redshift up 
to z ~ 1.1, with a net bias (z p ho — z gcn) of 0.008 on average. 

For only 1.4% of the events do we obtain a catastrophic 
error (|Az|/(l + z) > 0.15) on the estimate of the red- 
shift, and for 10% an error larger than 0.07. The rate of 
catastrophic failures is reduced to 0.4% for spectroscopi- 
cally confirmed SNe la. 

The method is illustrated with SNLS data but can be 
adapted to other supernova samples. 
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Fig. 9. Top: bias on redshift as a function of the fitted red- 
shift with the additional SALT2 fit for three bins of intrinsic 
(generated) color c gcn : black squares for |c gen — m e \ < 0.1, 
filled red circles for c, 
for 1 



gen 



c > 0.1 and open blue circles 
< -0.1. The curves illustrate the 0.014(1 + z) 2 
statistical uncertainty. Bottom: bias on color as a function 
of the fitted redshift. Same color code as above. 



-gen 



7. Conclusions 

The study presented in this paper shows that SN la photo- 
metric redshifts can be obtained with an accuracy exceed- 
ing that of current host galaxy photometric redshifts. The 
method was applied to simulated light curves as well as to 
a sample of SNe la from the first three years of SNLS data 
for which the redshift was known from spectroscopy. The 
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